Toroidal bubbles with circulation in ideal hydrodynamics. A variational approach 
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Incompressible, inviscid, irrotational, and unsteady flows with circulation F around a distorted 
toroidal bubble are considered. A general variational principle that determines the evolution of the 
bubble shape is formulated. For a two-dimensional (2D) cavity with a constant area A, exact pseudo- 
differential equations of motion are derived, based on variables that determine a conformal mapping 
of the unit circle exterior into the region occupied by the fluid. A closed expression for the Hamilto- 
nian of the 2D system in terms of canonical variables is obtained. Stability of a stationary drifting 
2D hollow vortex is demonstrated, when the circulation is relatively large, gA 3 ^ 2 /V 2 <C 1. For a 
circulation-dominated regime of three-dimensional flows a simplified Lagrangian is suggested, inas- 
much as the bubble shape is well described by the center-line R(£, t) and by an approximately circular 
cross-section with relatively small area, A(£,i) <C {S |R'|d£) 2 . In particular, a finite-dimensional 
dynamical system is derived and approximately solved for a vertically moving axisymmetric vortex 
ring bubble with a compressed gas inside. 

PACS numbers: 47.55.Dz, 47.15. Hg, 47.10.+g, 47.32. Cc 



I. INTRODUCTION 



Vortex ring bubbles in water are like usual vortex rings 
with circulation, but the core is filled with air, thus they 
are also termed "air-core vortex rings" . The higher ve- 
locity fluid surrounding the core of the ring is at lower 
pressure than the fluid farther away due to the Bernoulli 
effect. Vortex ring bubbles can be generated in various 
ways naturally or artificially, and they are interesting ob- 
jects both from experimental and theoretical points of 
view. Amazing examples of the natural beauty are vortex 
bubbles blown by dolphins for amusement. Also whales 
sometimes blow ring bubbles that can reach several me- 
ters in diameter. In laboratory conditions, toroidal bub- 
bles can be created relatively easy by an air jet that is 
rapidly opened and closed at the bottom of a water tank, 
as in the early experiments by Walters and Davidson [l| . 
The toroidal bubbles with circulation were formed as the 
result of gravity- induced topological transformation of an 
initial large spherical bubble, when a "tongue" of liquid 
penetrated the bubble from below. This appears to be 
a generic way of the creation of bubble rings (see, e.g., 
[Ill, and references therein). Such spherical bubbles 
may also be produced in nature, for instance, by un- 
derwater explosions. When formed the toroidal bubbles 
propagate upwards with an increasing diameter. Amus- 
ing examples of the generation and dynamics of vortex 
ring bubbles - or "silver rings" - may be found at the web 
site 0. 

The first attempts to describe the dynamics of the vor- 
tex ring bubbles analytically have been made a long time 
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ago (see 0- 0- El • and references therein) . It is clear that 
the most general and realistic theoretical consideration 
should be based on the Navicr-Stokcs equations, and thus 
is a very complicated nonlinear free-boundary problem in 
three-dimensional (3D) space. However, in many cases 
the inviscid approximation, based on the Euler equa- 
tions, may provide useful results. Inviscid flows belong 
to the class of conservative dynamical systems and thus 
are more easily studied by Hamiltonian and Lag rangian 

methods 0,SH El El El El El 13 El m EI Mm- 

With these methods, it is possible to simplify the anal- 
ysis considerably and make it more compact, especially 
for irrotational flows, when the original 3D problem re- 
duces to effectively a 2D problem on the free surface 

HE! El El El El El Ella 

The vortex ring bubble is a special example of a general 
class of fluid dynamical problems involving the free sur- 
face separating the fluid and air (or generally two differ- 
ent fluids). The aim in this paper is to develop a Hamil- 
tonian formalism for these systems, and the symmetric 
vortex ring bubble will be treated as a particular exam- 
ple. In other terms, we consider the question about the 
principle of least action for a general toroidal bubble. The 
corresponding Lagrangian functional is shown to possess, 
besides quadratic (inertial) terms on generalized veloci- 
ties, gyroscopic terms (of the first order on generalized 
velocities) . The gyroscopic terms are proportional to the 
constant circulation along linked contours. This prop- 
erty makes the toroidal bubble similar to a vortex fila- 
ment, if the circulation is large. We should emphasize 
that our approach, which is based on inviscid flows, can 
naturally not describe the topological transformation of, 
e.g., a raising spherical bubble into a vortex ring bubble 
as briefly discussed above. 

Having obtained a general variational formulation, we 
derive various approximations with reduced number of 
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degrees of freedom. First of all, we consider the ex- 
act reduction corresponding to purely 2D flows around 
a cavity. In this case it is possible to express the La- 
gra ngian in terms of the so-called conformal variables 
El El 13 El 13 HJ. For 3D flows, we do not have 
an exact explicit expression for the Lagrangian, but ap- 
proximations are possible. Such approximate dynamical 
systems take into account only the most relevant degrees 
of freedom of the original system. In this way we have 
obtained an approximate Lagrangian for a relatively long 
and thin toroidal bubble. For an axially symmetric ris- 
ing and spreading vortex ring bubble our variational ap- 
proach provides a finite-dimensional approximate system 
that is a generalization of the model discussed by Lund- 
gren and Mansour 

The remaining of the paper is organized as follows: 
In Sec. II we derive the general variational principle 
for the bubbles. We consider as an illustration the ex- 
ample of bubbles without circulation, before we derive 
the general Lagrangian for bubbles with circulation. The 
two-dimensional hollow vortex is considered in Sec. Ill, 
while the three-dimensional hollow vortex with a general 
toroidal shape is considered in Sec. IV. As a specific ex- 
ample we investigate the dynamics of the axisymmetric 
vortex ring bubble. Finally, Sec. V contains our sum- 
mary. 



II. VARIATIONAL PRINCIPLE FOR BUBBLES 

A. Hamiltonian structure of equations of motion 

It is well known that a class of irrotational solutions 
exists in the framework of ideal hydrodynamics. Such 
solutions describe potential flows with zero curl of the 
velocity field at any moment of time in the bulk of the 
moving fluid. If the liquid is also incompressible (with 
the unit density, for simplicity) then the investigation of 
non-stationary irrotational flows in a space region D with 
the free surface £ can be reduced to the consideration of 
Hamiltonian dynamics of the surface j^,|3QI3EH3 
L3I3I3] In this formulation, the shape of the surface E 
itself and the boundary value \& of the velocity potential 
are the dynamical variables determining the state of the 
system. The velocity potential of incompressible fluid 
satisfies the Laplace equation in the bulk of the fluid 



V(r, t) = V$, A$ = 0, $| s = 



(1) 



Besides the free surface E, in the general case the to- 
tal boundary of the region D has other pieces, which 
consist of infinitely far points of D and/or some walls 
(such as surfaces of submerged bodies). For the remain- 
ing boundary conditions for <I>(r), we will suppose that 
$(r) vanishes at infinitely far points, and employ the 
no-penetration condition at a motionless wall W (in par- 
ticular, it can be the bottom with arbitrary profile), if it 
is present: 



where N is a normal vector on the wall. 

The equations of motion for £ and W take the following 
form: 



£ = V n = (n-V$)| S) 

• <9$ V 2 
-* = I -E- — + — +.gz + P(V) 
on 2 



-V* + ^+gz + P(V) 



(3) 



(4) 



This set of equations describe, e.g., the dynamics of a 
bubble or a void submerged into a fluid. £ is the speed of 
surface motion along the normal unit vector n directed 
inside the bubble, V n is the normal component of the 
velocity field, and \E' is the total time derivative of the 
boundary value of the potential \& on the moving surface 
d [w| . The (normalized to the fluid density) pressure 
P(V) of the gas inside the bubble (the gas is considered 
as approximately massless and adiabatic) depends on the 
total volume of the bubble, 



V = 



-iy"(r-n)dS, 



(5) 



where dS is an element of the surface area. The vertical 
Cartesian coordinate z is measured from the horizontal 
plane where the pressure is zero at the equilibrium. The 
gravitational acceleration is —ge z . Thus, at the horizon- 
tal surface of the fluid at the atmospheric pressure (for 
instance, at the see level), z — z* w -10 m. Eq.@ is 
simply the kinematic condition, and Eq.Q follows from 
the Bernoulli equation for non-stationary potential flows 

m. 

It is possible to verify that the right-hand sides of the 
equations ©-© have the form of variational derivatives 



JH{E,*} 

££ ' 



(6) 



where the Hamiltonian 7i{£, ty} is the sum of the kinetic 
energy of the moving fluid, the internal energy of the 
compressed gas, and the potential energy in the uniform 
gravitational field (all the quantities are normalized to 
the fluid density), 

W = I [ (V$) 2 dr + £{V) + 9 - [ (e z ■ n)z 2 dS. (7) 

2 J D 1 JY, 

Here the adiabatic relation between the internal energy 
and the pressure is used, 



S'(V) = -P(V). 



(8) 



The derivation of the equality STC/S^ = V n is easy. In- 
deed, due to the equations Q - © one can write, 



SH\ 



<5£=0 



D 



V n S^ dS. 



(W-N)|vk = 0, $| 



0. 



(2) 



The calculation of 5H/SH is a bit more involved. It con- 
sist of two parts. First, due to the variation <5E (in the 
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normal direction) of the integration domain D without 
changing the potential $ inside, the following terms arise 



6J& 



5*=0 



2 



gz + P(V) <5£ dS. 



The second part comes from the condition that the value 
\& on the new (S + <5E)-surface must remain the same as 
on the old E-surface. To satisfy this requirement, the po- 
tential $ on the old boundary should acquire the change 
— (<9$/9n)5E. Therefore the second term is 

*« (a W = / S • «s = - I sn dS 



Sty 



dn 



The comparison of the sum of two these parts with the 
Eq-Q gives the second equation from Eq.JfjJ). 

Finally, we note that surface tension can be accounted 
for by simply adding the term a Jy,dS, which is the sur- 
face energy, to the Hamiltonian (Jj. Here a designates 
the surface tension coefficient (divided by the fluid den- 
sity). 



B. Variational principle 

1. Bubbles without circulation 

It is observed that in the simplest case, when the po- 
tential $ is a single-valued function, the equations of 
motion JSJ follow from the variational principle SA = 
5 J Cdt = with the Lagrangian 



£ = [ ty± dS-H{Z,ty}. 



(9) 



This expression is written in the invariant form that does 
not depend on the choice of the parameterization for the 
surface shape. Practically, this choice is dictated by ge- 
ometry of a given problem. For instance, the parameter- 
ization z = r)(x,y,t) is commonly used to study waves 
on the sea surface. It is clear that due to the equality 
ty~E dS — tyrjt dxdy the functions r)(x, y) and ty(x , y) form 
a pair of canonically conjugated variables But if 

we want to study oscillations of a spherical bubble, the 
spherical coordinates r, 9, tp are more convenient. In this 
case the functions — ty(9, ip) and Q(9, ip) = r 3 (9, <p)/3 can 
be taken as canonical variables. 

As an illustrative example, we consider the case cor- 
responding to spherically symmetric flows with g = 
and with a constant external pressure P ex t- In this case 
the dynamical variables depend only on t, and we have 
the completely solvable conservative system for ty(t) and 
V(i), represented by the Lagrangian: 



C 



,ty 



sph. 



_ aV 1 / 3 — - £(V) - P ext V - 6V 2 / 3 , (10) 



where a = 3 1 / 3 (47r) 2 / 3 , and b = 3 2 / 3 (47r) 1 / 3 o- accounts for 
the surface tension. The equations of motion - the Euler- 
Lagrange equations - corresponding to this Lagrangian 



-V - aV 1/3 <F = 0, 



*-™* 2 -£'(V)-P e: 2 



6V 2 / 3 



3 V 1 / 3 



0. 



From here one can eliminate *F and obtain the equation 
of the second order for V: 

V 1 V 2 2b 

P ext + £'(V) + -^ m = 0. (11) 



aV 1 / 3 6 aV 4 / 3 



It is easy to show that Eq. l|ll|) is equivalent to the sim- 
plest variant of the Rayleigh-Plesset equation for the 
radius of a spherical bubble (see [2^ and references 
therein), using V = (47r/3)i? 3 , 

RR+lk 2 = -P ext -S'(V)-^. 
z H 

Since the Lagrangian (|1(J|) does not depend explicitly on 
time, the system possesses the energy integral 

a yi/3_ + £(y) + PextV + 5v 2/3 = E ^ 
Therefore the solution of the Ea. Hllfl is determined by 

r v dv 

t = 



v » ^2aV!/3[£; _ £ (V) - P ext V - bW 3 } 

where E$ and Vo are arbitrary constants. If £(V) = 
(the bubble may contain no gas), then the above expres- 
sion describes a spherical cavity collapse for P ext > 0, as 
well as possible cavity formation for negative P ex t- 

More complex spherical bubble dynamics with a time- 
dependent P ex t(t) is governed by the Lagrangian ltTU|) 
as well, however, we do not have analytical solutions 
for that case. For instance, the dependence P ex t(t) = 
Po + P s cos(uit) is related to the problem of single bub- 
ble sonoluminescence |22t . where u> is the frequency of a 
(relatively long) standing ultrasound wave. 



2. Toroidal bubbles with circulation 

The variational formulation becomes more complicated 
in the case when the free surface E : £) i— ► R = 
(X,Y,Z), with < i? < 2tt and < ^ < 2tt, is topo- 
logically equivalent to a torus, and the circulation of the 
velocity along linked contours takes a nontrivial value T. 
Now the potential $ is a multi-valued function, 



$ = cf) + (r/27r)0, 



(12) 



where <f> is the harmonic potential determined by a single- 
valued boundary function -0(t9,£), and the velocity field 
created by the multi- valued harmonic function 9 has zero 
normal component on the free surface. The important 
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point is that the potential 9 is completely determined 
by the shape of the toroidal bubble. The multi-valued 
boundary function 0($, £) associated with the potential 9 
increases by the value 2tt as the coordinate i? acquires the 
increase 27r. The kinetic energy of the flow is represented 
as the sum of circulation-induced energy and the energy 
associated with the motion of the bubble. In the general 
form, we have the following expression: 

1 



It is interesting to note that for circulation-dominated 
configurations (it is important that the gradient of W 
along the surface should not be equal to zero at any 
point of S), a similarity exists between a hollow vortex 
tube and an ordinary toroidal vortex sheet. Indeed, the 
dynamics of a toroidal vortex sheet in a fluid without 
free boundary is governed by the Lagrangian (see, for 
instance, |23| and references therein) 



K = r 2 /C c {£} + - J J gx(s 1 ,S2W(s 1 )ip(s 2 )dSidS 2 , £ vsh = IJ ([R c xR t ].R)d^(-K ¥ . sll .{R(v,0}, (18) 



(13) 

where s% £ S, S2 £ S, and Gs(si,S2) is a symmetric 
function completely determined by a given shape of the 
bubble. 

In order to have correct equations of motion for 
R($, £,<) and ?/>($,£,£) (the equations must be equiva- 
lent to Eq.©), it is necessary to include into the action 
A = J Cdt a term that will give the same contribution as 
the following term, 



where the vector function R(^, £) describes the shapes 
of individual vortex lines enumerated by the label v € 
[0,r], and is directed along the vortex lines on the 
sheet. On the other hand, when considering the hollow 
vortex tube with a strong circulation, we could use the 
function ^ as a coordinate on the bubble surface, instead 
of the coordinate and in that case the Lagrangian of 
the hollow vortex tube would take the alternative form: 



2tt 



dt / esds 1 = 



2vr 



(19) 



It is clear that this expression should be transformed by 
some integration in parts to a form where is not em- 
ployed, but only the derivatives ©t, 0£, and 0$ that are 
single-valued functions. As a result, we obtain that the 
Lagrangian for a hollow vortex tube can be written as 
follows 

£ = [ ip± dS-H{Z,ip} 
is 
L 



3 • 2tt 



R 



{[R ? x Ri]0^ 



+[R t x R fl ]0 5 - [R 6 x R^]0 t | M df. (14) 

Now we may identify the function with the coordinate 
1? and thus the two last terms are equal to zero. Also it 
is possible in general to express the potential ip as 



ip(s) = / M s (s, s)±(s)d§, 



(15) 



where the "matrix" Ais is the inverse of the "matrix" 
Qs, and thus exclude ip from the Lagrangian. Then we 
will obtain the Lagrangian of the form 

£ = ^ / (R • [R« x R t ]) dd di - n{£} 



+ 2j J M ^ s ^( s i) E ( s 2)dSidS 2 , (16) 

where the effective potential energy n{£} is the sum of 
the circulation-induced energy, the internal energy of the 
compressed gas inside the bubble, the gravitational en- 
ergy of the bubble, and the surface energy, 

II{£} = r 2 /C c {£}+£(V) + ^ / (e 2 -n)z 2 dS + a [ dS. 

(17) 



where H{R(^ , £)} is the total energy of the toroidal 
bubble. Thus, the only difference between Eci. (|18|) and 
Eci. (|19|) is in the Hamiltonians H v . s h. and Tt. In the limit 
of a "thin vortex tube" the Hamiltonians are almost iden- 
tical, inasmuch as the main contribution is due to the 
logarithmically large circulation-induced kinetic energy. 

In the general case a free surface may consist of sev- 
eral separated manifolds with nontrivial topology. All of 
these must be included into the Lagrangian in a similar 
manner. 



III. 2D HOLLOW VORTEX 

As application of the theory described in the previous 
section, let us first consider a 2D irrotational flow in the 
j/z-plane, with the circulation T = 2n"/ around a cavity 
having a finite area A — nr 2 . The 2D geometry allows 
us to employ the theory of conformal mappings to derive 
exact equations of motion for such a system. Conformal 
variables have been extensively used during recent years 
for analytical studies of waves on water surface, and for 
numerical simulations (see, for instance 0, 0, 0, 0, 
El HI). The syste m considered in this section has a 
set of additional properties in comparison with the usual 
surface waves. The presence of the circulation makes it 
similar to a vortex. At the same time, the hollow vortex 
possesses inertial properties and a potential energy in the 
gravitational field. For small values of the parameter fi — 
gr 3 /j 2 a stationary horizontal drift of the hollow vortex 
is possible with the velocity Va ~ ~gr 2 /(2^y) and with 
the shape close to circular. This motion is stable, as will 
be discussed below. Therefore the content of this section 
will serve as a basis for further simplified descriptions of 
3D circulation-dominated flows. 
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A. Conformal mapping 

We consider an infinite two-dimensional region D, 
which is topologically equivalent to the exterior of the 
unit circle. Our purpose is to obtain an expression for 
the kinetic energy of the irrotational flow with the circu- 
lation 2-777 around the cavity in the case of an arbitrary 
given shape of the surface and arbitrary given bound- 
ary potential \& (with the only condition W i— ► if? + 2tt^ 
after one turn along the boundary. Strictly speaking, 
this energy is infinite because of the divergency of the 
corresponding integral at the infinity. But this is not 
important for the equations of motion, inasmuch as the 
presence of an infinite constant term in the Hamiltonian 
in no way influences the dynamics. Therefore only the 
excess of the energy in comparison with some basic state 
is needed. As the basic state, we shall take the perfect 
circular shape of the boundary, with the radius r and 
purely azimuthal velocity field, inversely proportional to 
the distance from the central point. 

Since the velocity potential, $>(y,z,t), satisfies the 
Laplace equation $ ya + $ Z z = 0, which is conformally in- 
variant, it is natural to re- formulate the problem in terms 
of the conformal mapping of the unit circle exterior into 
the region D. This mapping is determined by an ana- 
lytical function C*( w t t) — y{w, t) + iz(w, t) of a complex 
variable w. The function C*(ti>, t) has no singularities at 
\w\ > 1 and behaves proportionally to w when w — > oo. 
Therefore the expansion of this function in powers of w 
contains no positive powers higher than 1. The shape of 
the free surface is given parametrically by the expression 

Y(ti,t)+iZ({>,t)=(4w,t)\ w=ei , 
o 

= C(<M)=Ci(iK*+ Yl trn(t)e im ». (20) 



m— — oo 



In the potential 'J we now explicitly separate the term 
7$, which is responsible for the circulation, 



V(<d,t) = 71? + 



+ OO 



im/i) 



Vv, 



(21) 



(22) 



The term 71? corresponds to the multi-valued harmonic 
function $o(^) = Re(— ^Lnw) with zero normal com- 
ponent of the velocity at the free surface. The single- 
valued function tf> is related to evolution of the boundary 
shape. It can be understood as the potential of surface 
waves. The excess of energy is the sum of two parts. The 
first part is due to the kinetic energy of the surface waves, 



E s 



= 2tt 



E 



|m||V>n 



m— — oo 



The other part arises in the circulational energy as the 
result of changing of the effective cavity size, and it is 



completely determined by the coefficient (1 : 
2tt7 2 



E-y ' 



111 



Now we have to introduce some necessary linear opera- 
tors El El El El El El to deal with boundary values 
of analytical functions. In Fourier-representation these 
operators are diagonal, 



H 



i sign(m), M m = \m\, 



= i(lTsign(m)). 

(23) 

Here the operator H is the Hilbert transformation. The 
operator pv _ ) excludes the Fourier-harmonics with posi- 
tive m, while excludes the harmonics with negative 
m. The following equalities will be used in the further 
exposition, 

M = -Hd&, P (T) = ~(l±iH), pW + pH = i. 

We have now prepared all the necessary tools, and we 
are able to write down the Lagrangian for a 2D hollow 
vortex in the conformal variables: 



4/ 



2tt7 2 



In 



C(0)e 



-i-d 



dO 

2^ 



* 2i 

(cc-cco 



7 

2r 2 



— j ipMt/} d&-- 



c-cy (c + co 

2i 



Ai 



dd 



+ f (\p(+\(e- 2M ) + Ap(-)(Ce 21 ")) d». 



(25) 



Here C, = dt(, £' = 9^C- Besides the obvious terms that 
were already explained in the previous discussion, in the 
Lagrangian £ C onf. there is the term proportional to the 
constant area of the cavity. Its presence provides mini- 
mum of the circulational part of the Hamiltonian on the 
perfect shape £ = re 1 ® + (o- To be punctual, we have also 
included the terms with the Lagrangian multipliers A and 
A in order to specify explicitly the analytical properties 
of the function £(#). 

The variation of the action with the Lagrangian l|25|) 
gives (after some additional transformations, see Ap- 
pendix A) the equations of motion for £($, t) and ^(1?, t), 



2iMi/;\ 

(py) 2 - (7 + ^') 2 



(26) 



2|C'| 2 

(C-C) 7 2 

2i 2r 2 ' 



+ (7 + 4>')H 



Hip' 

w 



(27) 



Of course, these equations can also be obtained directly 
by simply presenting the kinematic condition and the 
Bernoulli equation in conformal variables. 
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B. Canonical variables 

The Lagrangian l|25|) is written in terms of variables 
that are not canonically conjugated. For general pur- 
poses, e.g., of convenience for a nonlinear analysis, a pair 
of canonical variables can be found. As the canonical 
coordinate, we take the real function q{§, t) such that 

C(iM)=0(iM)e w , t) = (l + (28) 

After substitution into the Lagrangian (|25p. one can im- 
mediately obtain the expression for the corresponding 
canonical momentum t), 



.•7 



P 



= -yHq - p(~3 [</>(/3 - 0)\ - [^(/3 + 0)] (29) 



Now it is necessary to solve this equation for the potential 
ip in order to express the Hamiltonian in terms of q and 
p. The result of the calculations is (see Appendix B): 



-ip{q,p} = 



(p-jHq)(q-Mq) + H[(p-jHq)H(q- Mq)} 



(q - Mq) 2 + (H(q - Mq)) 2 



Thus, the Hamiltonian for a 2D hollow vortex is: 
H{q,p} = l/j{q,p}M^{q,p}dd 

2?T7 2 ( fq \ , <7o 



(30) 



+ ^2 / q{l-M)qdd 



+~ J (q sin + Hq cos??) 2 

x [(q' - Hq) cos $-{q + Hq') sin <&] d&, (31) 

where ^{(^p} should be taken from Ea. H30fl . and qo is 
the 0-th Fourier- harmonic of the function g($), 



qo = / q($) 



C. Linearized equations for small deviations from a 
circular shape 

Let us consider an initial stage of evolution of a hol- 
low vortex, starting from a nearly perfect - circular - 
shape (with the mean radius r), in the case of a small 
initial potential ip. It should be emphasized that the cir- 
cular shape is not a stationary solution in the presence 
of gravity. Therefore we are studying dynamics that can 
be, generally speaking, far from equilibrium. However, 
at least at sufficiently small times we will have 

C = re M b{d, t), 6(0, t) = b Q + Y, b- m e-^, b Q a 1, 

m,—l 

(32) 

and \b- m \ <C 1 if m > 1. From Eq.JUjJ) and Ea.lfSTjl we 

obtain the linearized (nonhomogeneous!) system 



i/>_ m = 2i — mip- m +—-2 (l-m)b- m -— (&_( m+ i)-<5 m i) 

where 5 m \ is the Kronecker <5. Eliminating ip, we have 
the set of equations 

2 

b-m = 2i—mb- m -\ — j (m -m)b- m +— m(6_( m+ i)-<5 m i) 

(35) 

A particular solution of this nonhomogeneous system is 



b-i = -f^ (exp(2ijt/r 2 ) - 1 - 2i-/t/r 2 ) , (36) 

(37) 



4i7 2 

b- m = 0, m > 1 



Applicability of the linearized equations implies that the 
velocity of the vortex motion is small in comparison with 
the velocity of rotation j/r. Thus, the parameter fi = 
gr 3 /"f 2 should be small, or at least the time t should be 
small, 

H\ sin(7i/r 2 )| < 1. 

If /i <C 1, then Eas. (|36[l - (|37|> are approximately valid for 
arbitrary t and describe the horizontal drift of the vortex 
with the mean velocity Vd = ~gr 2 /(2'y). The shape of 
the vortex in this limit remains almost circular. 

The general solution of the corresponding homoge- 
neous system is the linear combination 

b-m(t) = J2 G A eX P (~ U J ^) . ( 38 ) 
v 

where v is a discrete parameter. The dimensionless fre- 
quencies and the corresponding modes b^] n should 
be determined from the algebraic system 

( m -(n^ + m ) 2 )b^l = ^ m b% i+iy (39) 



f b- m — — 2rmp. 



(33) 



It is clear that in the case g ^ the modes b^L m are de- 
localized both in $-space and in m-representation. They 
can be partially classified by the number n of the last 
non-zero Fourier-harmonics. Therefore v = (n, ...), and 

(n - (f>("<-> + n) 2 ) b { l l ;r } =0. (40) 

From here we have v = (n, ±), and 

n(™>±) = -n ± y/K. (41) 

We see that regardless of the vortex size r and the 
value of 7, all the frequencies are real. From the other 
hand, we naturally expect an instability for sufficiently 
large r and/or small 7. But there is no contradiction 
at this point because for a large size and/or a small 
circulation the behavior of a coherent superposition of 
TO-delocalized modes with real frequencies is effectively 
exponential at small t. Therefore the linearized equa- 
tions for small deviations (from the circular, not from an 
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unknown stationary shape!) become un-valid very soon 
and the nonlinearity begins to play an essential role. It 
should be emphasized that in this subsection we have 
not considered an exact stationary configuration of the 
hollow vortex, since we do not have an analytical expres- 
sion for such solution. One can expect that a stationary 
shape strongly deviates from circular when the parameter 
\x increases considerably, and finally, above some critical 
value jU* ~ 1, a stable stationary solution does not exist. 
However, at small /i stable stationary solutions do exist, 
and the stationary shape is almost circular, inasmuch as 
the linearized equations for small deviations remains ap- 
proximately valid for arbitrary time and their solutions 
are stable. The questions of the determination of the 
critical value /i* and how instability develops are left for 
future investigations. 



IV. 3D HOLLOW VORTEX TUBE 

A. Simplified Lagrangian 

We proceed to a simplified consideration of a 3D thin 
and long closed hollow vortex tube with a smooth center- 
line 

R(Z,t) = {X(Z,t),Y(t,t),Z(Z,t)) 

and with approximately circular cross-section having a 
relatively small area A(£,t) = na 2 (£,t) < A 2 , where a is 
the radius of the cross-section, and A is the total length 
of the center-line, 



where 



A= <j)\R'\d£_. 



(42) 



This description should be good in most cases, since 
for large enough circulation (r 2 ^> gA 3 / 2 ) the local 
quasi-2D dynamics is stable with approximately circu- 
lar cross-section, and also for a straight thin 3D tube 
(r 2 > 47r 2 a (A J V) 1 ' 12 , see 5|) it is easy to demonstrate 
stability of longitudinal sausage-like perturbations. 

Assuming slow variation of A(£,t) along the curve 
R(£,i) and neglecting small distortions of the shape of 
the circular cross-section, we can give an explicit form 
of all terms in Eas. (|16() - (|17|l . As the result, a simplified 
Lagrangian with logarithmic accuracy can be written as 
follows, 



C=±f\n(¥.){A ± ) 2 \Rt\dt- 



\Ri 



-A\R 6 \d* 



~ f(R ■ [R e x R,]) d£,-^-j\n (^) iH^ldC 

s( <j^A\R e \d^j + g <j> ZA\Rf:\d£, 

-2^' 2 a j A^ 2 \R 6 \d^, (43) 



A ± = ^-^(R t -R c )/|R c | 2 , (44) 
R ± = R t R^Rt -R C )/|R ? | 2 . (45) 

This Lagrangian includes the most principal inertial 
effects that correspond to the dynamics of R(£, t) and 
A(£, t) on scales of the order of A. The interplay between 
the second order time-derivative (inertial) terms and the 
circulation-originated first order terms will result in os- 
cillations that are relatively fast if T is large. However, 
in the circulation-dominated regime this system has in- 
teresting solutions with these oscillations almost not ex- 
cited. Approximately such non-oscillating solutions are 
determined by the Lagrangian without the inertial terms. 
That means we have to find a minimum of the effec- 
tive potential energy n{R(£), A(£)} over A(£) with fixed 
R(£) and then substitute the minimum-providing config- 
uration A*(£) into n. The extremal configuration A*(£) 
is determined by the following coupled equations, 

+ P(V) + gZ{0 - Tr^aA^ 2 (£) = 0, (46) 



^A(£)|R'(£)K = v, 



(47) 



as observed from the Lagrangian l|43|) . At this point we 
meet a technical difficulty, since although Eq. I|4t)|) has the 
explicit solution 



.4 



-1/2 _ 



-in 



= ^ (V /2 CT + s/-KC 



^-{T 2 /2n)[P{V)+gZ(0]), 

(48) 

unfortunately Eq. (|47|l for V, with this expression for 
A*(0, is hard to solve exactly, except for the sim- 
plest case, P(V) = 0, when v4*(£) is not depending on 
the volume. Nevertheless, approximate methods may 
be used in many cases and corresponding approximate 
expressions for the effective Hamiltonian 7i*{R(£)} = 
n{R(£), ^4*(£)} of the bubble center-line can be obtained. 
The equation of motion will then have the general struc- 
ture as follows, 



[Re x Kt] 



1 STL 

r sr ' 



(49) 



which are similar to the equations of motion for a slen- 
der vortex filament in a fluid without bubbles jU I3| , 
however, with another Hamiltonian. 

For a hollow vortex tube without gas inside, when 
P(V) = 0, the dynamics of the center-line is described 
by the effective Hamiltonian 

Wf=°{R(0} = r^F(Z(0)|R'(£M, 
where the function F(Z) is defined as follows, 



F(Z) 



47T I 



Z 



C 



} 



c + vc 2 - z 

(50) 
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Here A ~ A<7 1//2 r _1 may be considered as approximately 
constant, and C 2 = 2ir 2 a 2 g~ 1 r~ 2 . The equation of mo- 
tion 149(1 for this case can be rewritten as 



R t = F'(Z)[e z x t] + F(Z)nb, 



(51) 



where t, b, and k are the unit tangent vector on the 
center-line, the binormal vector, and the curvature of 
the line, respectively. This equation is a generalization 
of the well known LI A (Localized Induction Approxima- 
tion) equation pi . 



B. Axisymmetric motion 

An obvious application of the Lagrangian (|43(l is for 
vertically rising and spreading axisymmetric vortex ring 
bubbles Q We need only three degrees of freedom to 
describe this dynamics, namely the vertical coordinate 
Z(t), the radius of the ring R(t) = (X 2 + Y 2 ) 1 ' 2 , and 
the total volume of the bubble V(t) = 2n 2 a 2 (t) R(t) . The 
corresponding finite-dimensional dynamical system is de- 
termined by the following Lagrangian, 



L 



ZRV 



(RV - RV) 2 
16n 2 R 3 



v 



Z 2 ) 



T 2 R 



i-\ ziu: u, n , — 

4 V lo s- V 
-S(V) + gZV - 2 3 / 2 naR 1 / 2 V 1/2 , 

where the constant coefficients 



(52) 



(A) = 
log. 



128tt j 



a 



(r) _ 



1287r 2 exp(-4), 



are used to improve the logarithmic accuracy. These co- 
efficients arise from asymptotic expansions of elliptic in- 
tegrals expressing the kinetic energy of a flow with a line 
source and with a vortex string on the ring center-line. 

Once the dependence P(V) is given explicitly, it is easy 
to write down and solve numerically the equations of mo- 
tion determined by the Lagrangian (|52|l . We used an 
approximate equation of state for the gas: 

P(V) =P (V /V) 1A . 

Examples of the solutions Z(t), R(t), and V(t), for sev- 
eral values of the circulation, are presented in Fig^ The 
solutions have an oscillating behavior with a drift. With 
such initial data, the amplitudes of oscillations for Z(t) 
and R(t) are smaller at larger T. However, oscillations 
for V(i) becomes stronger at highest T, that in practice 
should result in intensive sound irradiation (this natu- 
rally cannot be described by the present theory of in- 
compressible flows, but needs compressibility). 

The system <|52ll has 6D phase space, but its weakly 
oscillating solutions approximately correspond to the 2D 
(R, Z)-system that is obtained by neglecting the inertial 
terms (quadratic on the time derivatives) in Eq. (|52|l and 
minimizing the expression n(i?, Z, V) over V, 



n : 



^lnfc (r) ^- 

4 ln V Cl °s- V 



■ E (V) - gZV + 2 3 /' 2 iraR 1 / 2 V 1/2 . 




3 4 5 6 7 

time (s) 



FIG. 1: The evolution of the vertical position Z(t), the ra- 
dius R(t), and the volume V(t) of a vortex ring bubble, 
for different values of the circulation, V — 1.0, 1.5, 2.0, 
2.5, 3.0 m 2 /s. In these simulations the initial data were 
Vo = 0, Ro = 0, Z = 0, Vo = 1.0 m 3 , Ro = 2.0 m, and 
Z = -(r 2 R /(4V ) + P -2 1 ' 2 TTaR 1 /2 Vo 1/2 )/g, with g = 9.8 
m/s 2 , a = 7.5 ■ 10~ 5 m 3 /s 2 . The parameter P = 200 m 2 /s 2 
approximately corresponds to the initial pressure 2 atm. The 
curve Z(t) with the largest displacement, the curve R(t) with 
the weakest expansion, and the curve V(t) that reaches the 
largest values for V, correspond to the highest value of the 
circulation. 
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FIG. 2: Comparisons of the bubble ring trajectories in ZR- 
plane, corresponding to the approximations 1521 and 1541 . for 
different values of the circulation. Parameters as in Fig0 



Thus, we have to solve the equation dH(R, Z, V) jdV = 0, 

t 2 r 



4V 



P(V) - gZ + 2 1/2 7rcr-R 1/2 V 



? l/2,;-l/2 _ 



0, (53) 



and find from here an equilibrium value V*(P, Z). The 
slow dynamics is approximately described by the La- 
grangian 

L RZ = -IttTZRR - Il(R, Z, V*(R, Z)). (54) 

However, since Eq. (|5Hfl is in fact already solved for 
Z(R,V), it will be convenient to rewrite the Lagrangian 
(|54|l in terms of R and V: 



Lrv — 



2ttT /Y 2 R 



9 V 4V 

r 2 p, 

m 



P(V)-2 1 /2 7r(7j Rl/2 V -l/2j i? ^ 



(Cg^-) - £{V) - 2^aR 1 ' 2 V 1 / 2 



-V 



r 2 p 

4V 



+ P(V) - 2 1 / 2 naR 1 / 2 V- 1 / : 



(55) 



This dynamical problem is, of course, completely inte- 
grable. Phase trajectories in the (R, V)-plane are the 
level contours of the effective Hamiltonian 



Hrv — 



T 2 R 



4 v~ V~ log ' V 
+2 1/2 7rcjP 1/2 V 1/2 = E a = const. (56) 
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The comparison between the drifting solutions of the 
6D system (I52II and the corresponding solutions of the 
2D system (I54f) is given in Fig. [21 We observe that the 
agreement is very good in particular for the high values 
of the circulation. Furthermore, the overall evolution is 
qualitatively in agreement with the numerical results of 
Lundgren and Mansour 0, especially for large values of 



FIG. 3: Simulations for dimensionless quantities: g = 1, cr = 
0.025, T = 5, V = 4tt/3, R /a = 5, and Pn = 10, 100, 1000. 
The parameters Vn, Ro, Zq for the 6D approximation were 
calculated from the 2D approximation in order to minimize 
the oscillations. The difference between 2D and 6D approxi- 
mations is almost invisible for Z(t) and R(t). The trajectories 
with large Pn are in good agreement with Fig. 3 in Ref.0. 
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the parameter Pq, when the relative variation of the vol- 
ume is small (see Fig|3 where all the plotted quantities 
have been made dimensionless as in Ref.Q). 

It should be kept in mind that applicability of the ax- 
isymmetric (both 6D and 2D) models requires, besides 
2n 2 R 3 > V and T 2 > g{V/Rf /2 , at least one more con- 
dition: 



i?-V2 v i/2 < 



r 2 



23/2 7 



(57) 



This inequality makes constant-cross-area configurations 
be stable to longitudinal sausage-like perturbations [5j . If 
a vortex ring bubble violates this criterion in the course 
of motion, then it will be destabilized by surface tension 
and the development of the instability, together with vis- 
cous effects, will result in transformation of the ring into a 
closed chain of smaller bubbles connected by thin vortex 
filaments, as it is observed experimentally (concerning 
this phenomenon, Lundgren and Mansour |2| refer to a 
private communication by D. McSweeny (1990)). How- 
ever, we found with realistic dependences £ (V) that the 
Hamiltonian evolution of the spreading ring, described 
above, typically moves away from the instability crite- 
rion. Thus, to destroy the ring, viscous diffusion of the 
vorticity seems to be necessary 0, This will result 
in a slow decrease of the circulation along a contour just 
above the bubble surface. Ultimately, this circulation 
becomes too small to satisfy the stability condition (|57|l . 
the instability will then evolve and the ring will tend to 
break up. 



V. SUMMARY 

In this work we have developed a variational approach 
for the theoretical study of the ideal incompressible ir- 
rotational flows with a distorted toroidal bubble, in 
the case when the velocity potential is a multi-valued 
function. Using this method, we have derived exact 
pseudo-differential equations of motion for purely two- 
dimensional flows with circulation around a single cav- 
ity. Also we have suggested a simplified Lagrangian for 
a three-dimensional thin hollow vortex tube. As a sim- 
ple particular case, the axisymmetric vertical motion of 
a spreading vortex ring bubble with a compressed gas 
inside has been considered. Approximate solutions of a 
corresponding finite-dimensional dynamical system have 
been obtained. 



Acknowledgments 

These investigations were supported by INTAS (grant 
No. 00-00292). The work of V. R. also was supported by 
RFBR and by the Russian State Program of Support of 
the Leading Scientific Schools. 



APPENDIX A: VARIATION OF J £. conf dt 

Here we extract the equations of motion (|26|l and l|27() 
from the expression (|25|) . First of all, variation of the 
action J C con {,dt by Sip gives the kinematic condition in 
terms of the conformal variables 



(CC 7 - ceo 

2i 



Now a standard procedure is to divide this equation by 

IC'I 2 , 

C C _ 2iMip 

c^c 7 " W 

and apply the operator p(~) which excludes Fourier- 
harmonics with positive m. As the result, we get Ea. l|26() 

e = C 'p(-) ( 2i ™A 



IC'I 



The variation of J £ con fdt by and subsequent exlu- 
sion of the Lagrangian multiplier A gives the equation 



P<->((-C( 7 + « + (^+«-£k' 



7 le 



21 



2 § 



e ** =0. 



The statement p(~) / = means that / contains only 
harmonics with positive m. The function C'e does not 
contain harmonics with negative m. Therefore the equal- 
ity P(~> (ft, 1 e % ®) = is true. In our case this results in 



^(-cc^ + ^ + ^^ + ^-^ici 2 



7 2 i('e M 



2 Ci 



0.(A1) 



It is easy to check that 



7^ iCe^ \ _ 72 



Now we have to take the real part of Eq. I|A1|) . Using the 
property P^ + p(+) = 1, we get the equation which in 
fact is solved for ip, 

(p<-»{( 7 + ^CfP<-.(^)} + c.c.)4. 

After simplification that uses the equality 
H{2iP' Mt};) = (V>') 2 - (Mip) 2 , 
we finally obtain Eg. (1271) . 



11 



APPENDIX B: EXPRESSION (l30l) 

To solve for ip the equation 

p = jHq - P ( ") [ip{/3 - i/3')} - P (+) [ip{{3 + 0% 

we perform the following steps. First, let us separate the 
harmonics with m < 0: 

P (-) [ip(J3 - i/3') +p- jHq] =ia = 2iP ( ~ ) a, 

where a is an unknown real quantity constant on Now 
we use the same trick as in the Appendix A, i.e., we 
multiply the above equation by the function j3 + i(3' , 

pi-) tyip _ ipf + (p _ 7 # g )(^ + 0')] 
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